The geometric structure of the Landau bands 
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We have proposed a semiclassical explanation of the geometric structure of the spectrum for 
the two-dimensional Landau Hamiltonian with a two-periodic electric field without any additional 
assumptions on the potential. Applying an iterative averaging procedure we approximately, with 
any degree of accuracy, separate variables and describe a given Landau band as the spectrum of 
a Harper-like operator. The quantized Reeb graph for such an operator is used to obtain the 
following structure of the Landau band: localized states on the band wings and extended states 
near the middle of the band. Our approach also shows that different Landau bands have different 
geometric structure. 

PACS numbers: 73.43.-f, 03.65.Sq 



The standard theories of the integer quantum Hall ef- 
fect are modelled by the Hamiltonian for a charged 
particle in a heterostructure with a uniform magnetic 
field; they are based on the following assumption concern- 
ing the structure of the spectrum: the spectrum consists 
of broadened Landau levels (Landau bands) generated 
by extended states near the middle of each band and 
by localized states in the band wings; then pinning the 
Fermi level by the localized states causes the occurrence 
of plateaus in the conductance curve (see Q for details). 
This specific spectral structure is usually related to the 
presence of impurities or other defects in the quantum 
Hall system ||. Indeed, for the Gauss or Lloyd ensem- 
bles of impurities the density of states has the behavior 
required for the explanation of the quantum Hall effect 
||. Moreover, using ideas from || it can be rigorously 
shown that the desired spectral structure is displayed by 
the two-dimensional Landau Hamiltonian with substitu- 
tional point impurities arranged at the nodes of a lattice, 
assuming a quasi-periodic distribution of the coupling 
constants ||. The same holds for point impurities with 
random positions on the system plane and with random 
distribution of the coupling constants . 

The aim of the present article is to show that this spe- 
cial spectral structure can also be verified by semiclas- 
sical approximation for a purely periodic Hamiltonian 
with a uniform magnetic field without any additional as- 
sumptions on the potential. Note that a semiclassical 
approach to the quantum Hall effect has been consid- 
ered in |^] for a Hamiltonian with a disordered poten- 
tial, within the framework of percolation theory on the 
Chalker-Coddington network. The results of ||] confirm 
the intuitive picture which relates the extended states 
in the band center with the central drift of an open cy- 
clotron orbit, and the localized states in the band wings 
with closed orbits. This intuitive picture was rigorously 
justified in |10| using Maslov's canonical operator JL1|. 



fllOf which describe in detail the semiclassical spectrum 
of the periodic Landau Hamiltonian and show, in partic- 
ular, the desired spectral structure. 

We use the Landau gauge for the Hamiltonian of a 
magneto-Bloch electron, 

H = !-(Hdi + (eB/ch)x 2 ) 2 - d 2 2 ) + V(x u x 2 ), 

where the potential V is periodic with respect to the lat- 
tice A generated by the vectors ai = (L, 0) and a 2 ; we as- 
sume that V is a real-analytic function. Introducing the 
cyclotron frequency lo c = \eB\/cm, the magnetic length 
hi = (h/muic) 1 / 2 , and passing to dimensionless potential 
energy v = V/ max V and to coordinates X = x/L, we 
rewrite H in the form H — mL 2 ui 2 H° where 



(A + a 2 ) 2 + p 2 2 



+ s v v{X u X 2 ). 



Below we offer some consequences of the results from 



Here Pj = —iesd / 9Xj and the dimensionless quantities 
£b = {l/L) 2 and ey = Sb max |V|//iw c will be consid- 
ered as small parameters in our spectral problem. Let 
us estimate these parameters in a typical situation: for 
a field strength B ~ 10T the magnetic length is of order 
lOnm; the characteristic length L for periodic arrays of 
quantum dots or antidots in GaAs is of order 200-500nm, 
such that £b ~ 10~ 3 ; and using for GaAs the electron 
effective mass m — 0.067m e we estimate Hu) c as 15meV, 
implying that for V with values in the range 1-25 meV 
fl3f we have ey < £b- 

The main difficulty of the problem under study con- 
sists in the non-integrability of the corresponding clas- 
sical system. At first sight, the presence of the small 
parameter ey suggests the use of perturbation theory; 
but it is readily seen that this will only be possible if we 
impose additional relations between sb and ey, because 
eb is also small. Moreover, using perturbation theory we 
can obtain only rather rough spectral information. 
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The smallness of sb, however, will allow us to describe 
the fine structure of each Landau band by using a semi- 
classical approach. 

The trajectories of the classical system with the un- 
perturbed Hamiltonian ^ ((Pi + X 2 ) 2 + P| ) on the plane 
(Xi,X2) are the cyclotron orbits with radius y/2I and 
angle <p centered at a point (2/1,2/2)- Therefore, it is rea- 
sonable to introduce new canonical variables: the gener- 
alized momenta /, y\ (or p, yi) and the positions ip, y 2 
(or q, 2/2) according to the formulas 

Xi=q + Vi, Pi = -V2, X 2 =p + y 2 , P2 = -q 
p = V2I cos ip, q = V2J sin ip. 

In these variables the classical Hamiltonian H° takes the 
form H° = I + evv{y/2I\ sin ip + yi, y/2I cos ip + y 2 ). 

This representation suggests that we average Hq over 
the angle <p, and describe the central drift of the cyclotron 
orbits by means of the averaged Hamiltonian 

2-7T 

H m (I,yuV2;e v ) = ^J H°dp (J) 


= I + evMy/-2IAy)v( yi ,y2), 

where Jo denotes the Bessel function of order zero; here 
the second term on the right-hand side results from the 
action of the pseudodifferential operator Jo(\/ — 2JA y ) 
on v. Next, it can be shown that there is a canonical 
change of variables (p',q',y') = {p,q,y) + 0{ey) such 
that 

H°(p,q,y;ev) = (±( P ' 2 + q' 2 ),y' ; e v ) +0(4), 

but the error estimate 0{e v ) is not good enough to de- 
scribe the desired fine structure of the spectrum of H°. 
The crucial point, then, is the following assertion which 
is proved in pToj j by suitably iterating the averaging pro- 
cedure: 

There is a canonical change of variables (P, Q, Y) = 
(Pj^jY) + 0(sv) such that in the new coordinates 
(P, Q,Y) the classical Hamiltonian system, H , be- 
comes integrable modulo 0(sy). More precisely, 
H°(p,q,y;ev) = H a (^{P 2 + Q 2 ),Yi,Y 2 ;e v ) + 0(ey), 
where the first term is real- analytic and all the terms are 
two-periodic in Y. 

It is worthwhile to note that the first step of the men- 
tioned iteration procedure gives just the function H av . 

Now we quantize 7i° considering the pairs {P, Q} and 
{Yi, Y2} as pairs of canonically conjugated variables, i.e. 
we set P = —iesd/dQ, Y\ = —ieBd/dY 2 . By the Corre- 
spondence Principle, a canonical transformation in classi- 
cal mechanics should correspond to a unitary transforma- 
tion of the quantum mechanical state space. In our case 
this expectation is expressed in the fact that the spectra 
of the operators H° and Ti. coincide semiclassically, up 



to 0(e^ + Sy ). Now, since P and Q commute with Yj 
(j = 1,2), the quantum Hamiltonian 7i° commutes with 
the Hamiltonian of the harmonic oscillator \(P 2 + Q 2 ). 
Hence the eigenfunctions \& of H° may be chosen in the 
form ^(Q,Y 2 ) — i/j n (Q)(p n (Y 2 ), where ip n are the eigen- 
functions of the harmonic oscillator corresponding to the 
eigenvalues E n = (n+l/2)eg, n — 0, 1, . . ., and the func- 
tions ip n satisfy the equation 

H n <p n = E<p n . (2) 

Here H. n is a pseudodifferential operator obtained as the 
quantization of the classical Hamiltonian TC n (Yi,Y 2 ) = 
H°(E n ,Yi,Y 2 ,ev)- As mL 2 uj 2 EB = huj c , the number E n 
is nothing but the n-th Landau level. Thus the spectrum 
of each operator TL n arises from broadening the n-th Lan- 
dau level into a band under the influence of the periodic 
potential. More precisely, in our approach each Landau 
band is described by the one-dimensional equation (Q), 
which depends on the band, and the original spectral 
problem is reduced to a family of one-dimensional spec- 
tral problems. This reduction makes it now possible to 
describe the fine structure of each Landau band using es- 
tablished methods of semiclassical approximation, since 
we only have to deal with an integrable Hamiltonian. 

In view of the established periodicity, 7i n can be 
treated as a Hamiltonian on the two-torus; therefore, 
only the following three types of trajectories occur for 
the corresponding Hamiltonian system: (1) contractible 
closed curves and extremum points, (2) non-contractible 
closed curves, (3) separatrices and saddle points. The 
concept of the Reeb graph flij has proved a very useful 
tool in describing the topological structure of the space 
of trajectories. Namely, points of the Reeb graph are in 
a one-to-one correspondence with connected components 
of the level sets of the Hamiltonian, in such a way that 
separatrices correspond to the branching points of the 
graph. If one represents H n as the height function of a 
suitably deformed torus, then construction of the Reeb 
graph becomes especially evident, which is illustrated in 
Fig. and[l]B. . 

Obviously, contractible trajectories on the torus are 
covered by families of closed trajectories on the plane 
(Yi,Y 2 ), and non-contractible ones are covered by fam- 
ilies of non-closed (but periodic) trajectories on this 
plane. Hence, to obtain the semiclassical approximation 
to the spectrum of 7i n , i. e., to obtain the semiclassical 
structure of the n-th Landau band, one has to quantize 
only the closed trajectories, using the Bohr-Sommerfeld 
rule; thus we obtain the so-called quantized Reeb graph 
(Fig. |l|C). Returning to the original coordinates (P,X) 
we see that the closed trajectories of 7i n generate two- 
dimensional tori. It is well known that the semiclassical 
eigenfunctions corresponding to these tori decrease expo- 
nentially outside the classically allowed regions (i. e. the 
projections of the tori onto the X-plane), hence corre- 
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FIG. 1: Two examples of the geometric structure of the Landau band: (A) Level curves of a Hamiltonian on the torus and 
their realization through the height function; (B) Trajectories on the torus: (1) contractible closed ones, (2) non-contractible 
closed ones, (3) separatrices, and the Reeb graph of the Hamiltonian; (C) Applying the Bohr-Sommerfeld quantization; (D) 
Relationship between the Reeb graph and the Landau band 



spond to localized states. On the other hand, open tra- 
jectories in the Y-plane correspond to two-dimensional 
cylinders in the phase space (P,X). Since those are 
non-compact the corresponding states must be extended; 
they can be constructed using Maslov's canonical opera- 
tor 0. 

Our considerations show that the nth Landau band in 
the spectrum of H is structurally completely determined 
by the Reeb graph of TL n (see Fig. |l]D). It is evident 
from this correspondence that, for a generic two-periodic 
potential V, the "highest" and the "lowest" points of the 
Reeb graph correspond to extremal points of H n . This 
means that '"boundary vertices" typically correspond to 
closed trajectories such that, consequently, the wings of 
the Landau band are filled by semi-classically localized 
states. On the other hand, the "middle region" of the 
Reeb graph contains loops or nodes which correspond to 
the open trajectories or separatrices such that the central 
part of the Landau band consists of extended states. 

Using the previous considerations one can describe the 
spectrum globally. To do this, one has to construct the 
Reeb graph of H°(I, Y\,Yz,Ey) for each /; as I runs 
through [0, +oo), the Reeb graph is transformed and 
traces out a certain surface (the Reeb surface of the po- 
tential v). We illustrate this by considering the potential 
v(X 1 ,X 2 ) = cosA"i + cosX 2 + 2cos(Xl + X 2 ). Using 



Eq. d), we get 

H av (I, Yi , Y 2 ) = / + e v ( Jo (V2I) (cos Y x + cos Y 2 ) 
+ 2J (V47) cos(Yi +Y 2 )). 

Use first 7i av instead of H°, then the Reeb surface has 
the form shown in Fig. |^. Quantizing I by the Bohr- 
Sommerfeld rule I = E n we select a discrete subset of 
the Reeb graphs, which can be obtained by cutting the 
Reeb surface by the plane I — E n . For different indices n 
we obtain different types of Reeb graphs, some of which 
are shown in Fig. 0. 

It is apparent that in all these cases the Landau band 
has the asserted structure. For some n, the domain of 
infinite motion may be very small and tend to a point 
(these points are emphasized in Fig. 2), but it always ex- 
ists. The inclusion of the correction TiP — H av may lead 
to the appearance of new vertices and new edges in the 
Reeb graph near the branching points; but these edges 
have length 0(e v ) and hence do not change the asserted 
spectral structure. Clearly, spectral regions, correspond- 
ing to different Landau bands may intersect. 

To summarize, we have proposed a semiclassical ex- 
planation of the geometric structure of the spectrum for 
the two-dimensional Landau Hamiltonian with a two- 
periodic electric field. Applying an iterative averaging 
procedure we approximately, with any degree of accu- 
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logically different graphs and, therefore, can have differ- 
ent structure. Such a distinction of the geometric struc- 
ture of the different Landau band was established numer- 
ically in jL5| for various kinds of potentials. 

In carrying out this work we had useful discus- 
sions with S. Albeverio, J. Avron, E. D. Bclokolos, 
V. S. Buslaev, P. Exner, M. V. Karasev, E. L. Ko- 
rotyaev, V. A. Margulis, A. I. Neishtadt, L. A. Pastur, 
M. A. Poteryakhin, A. I. Shafarevich, P. Yuditskii, and 
J. Zak. To all these people we express a gratitude. 
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FIG. 2: The Reeb surface for the potential v(Xi,X 2 ) = 
cos(Xi) + cos(X 2 ) +2cos(Xi + X 2 ) 




FIG. 3: The structure of different Landau bands for the po- 
tential v(Xi,X 2 ) = cos(Xi) + cos(X 2 ) + 2cos(Xi + X 2 ) 

racy, separate variables and reduce the original spectral 
problem to an infinite family of one-dimensional eigen- 
values problems for Harper-like operators; each of these 
operators is the quantization of a classical Hamiltonian 
having a torus as the phase space. Such a reduction 
provides a convenient tool for describing the geometric 
structure of each Landau band. Namely, the space of 
trajectories of each Hamiltonian on the torus is repre- 
sented by the so-called Reeb graph. The projection of the 
Bohr-Sommerfeld quantization of the graph onto the en- 
ergy axis gives the corresponding Landau band. Clearly, 
open edges of the graph represent contractible trajecto- 
ries; quantization of these trajectories leads to localized 
states in the Landau band. Therefore, the band wings 
are filled by the localized states. On the contrary, open 
trajectories are always represented by closed edges, which 
lie in the middle of the graph; they generate the extended 
states near the middle of the Landau band. 

It should be particularly emphasized that different 
Landau bands are described, generally speaking, by topo- 
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